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Abstract 

We consider a hydrodynamic model of swarming behavior derived from the ki- 
netic description of a particle system combining a noisy Cucker-Smale consensus 
force and self-propulsion. In the large self-propulsion force limit, we provide ev- 
idence of a phase transition from disordered to ordered motion which manifests 
itself as a change of type of the limit model (from hyperbolic to diffusive) at the 
crossing of a critical noise intensity. In the hyperbolic regime, the resulting model, 
referred to as the 'Self-Organized Hydrodynamics (SOH)', consists of a system of 
compressible Euler equations with a speed constraint. We show that the range of 
SOH models obtained by this limit is restricted. To waive this restriction, we com- 
pute the Navier-Stokes diffusive corrections to the hydrodynamic model. Adding 
these diffusive corrections, the limit of a large propulsion force yields unrestricted 
SOH models and offers an alternative to the derivation of the SOH using kinetic 
models with speed constraints. 
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1 Introduction 

There is a considerable literature devoted to the observation and understanding of systems 
of swarming agents. Examples of such systems in nature are fish schools (6j [39], bird 
flocks [21 HQ] , insect swarms [TH HE] or migrating cell assemblies [55] (see also the reviews 
[US El])- Some simple inanimate physical systems also exhibit collective behavior [291B7] . 
Many of the models proposed in the literature are 'Individual-Based Models (IBM)'. They 
consist in following the dynamics of the agents and their interactions over the course of 
time. The 'three zone model' of Aoki [H EQ] postulates that interactions obey long- 
range attraction, short-range repulsion and medium-range alignment. Vicsek et al. [53] 
have proposed a simplified version of this model where particles move at a constant 
speed and interact through alignment only. In spite of its simplicity, the Vicsek model 
exhibits complex features which have triggered a large literature (see [51] for a review and 
references). On the other hand, the Cucker-Smale model [22] is based on a large-scale 
velocity consensus formation and does not impose any constraint on the particle speed. 
The Cucker-Smale model has triggered considerable mathematical activity [T5l EU |35j [361 
E3 HH]. Many other kinds of IBM's of collective motion can be found and it is impossible 
to cite them all (see e.g. [TBI H3 EH1 H2] and the review [M])- Additionally, comparisons 
of models with data can be found e.g. in [H [5j [38]. 

IBM's are very successful but become computationally intensive for large systems. 
For this reason, macroscopic models of fluid type have been proposed in the literature. 
Macroscopic models of collective motion have been derived from heuristic rules and sym- 
metry considerations in [1H1 US]- The rigorous derivation of continuum models usually 
starts from a statistical version of the IBM, the so-called kinetic model. Kinetic models of 
collective motion have been derived in [TTJ [121 E5] for various versions of the Cucker-Smale 
and Vicsek models. The convergence of the kinetic Cucker-Smale model to the kinetic 
Vicsek model is shown in [10]. In [71 [8], a Boltzmann type kinetic model has been proposed 
for binary collision processes which mimics the Vicsek alignment dynamics. In the same 
spirit, a Boltzmann-Povzner type approach which mimics the Cucker-Smale process and 
its fluid limit has been developed in [32]. In [7J Ej, a hydrodynamic model for the binary 
version of the Vicsek interaction is derived from the kinetic model under an assumption 
of weak anisotropy of the velocity distribution function. In [15] a direct passage from 
the Vicsek IBM to a fluid model is attempted. The first derivation of hydrodynamic-like 
equations from the mean-field kinetic version of the Vicsek model has been performed in 
[26] . Further elaboration of the model can be found in [211 [23 [331 EI]- Diffusive correc- 
tions to the model of [2Ej have been derived in [2H] and bear analogies with the model 
proposed in [13]. Other kinds of macroscopic models can be found in [HI EH HH [501 El]- 

The aim of the present work is twofold. The first objective is to give evidence of a 
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phase transition from disordered to ordered motion in a hydrodynamic model of socially 
interacting agents with self-propulsion. Specifically, we want to emphasize the role of the 
self-propulsion in the emergence of the phase transition. Such evidence has been given for 
the first time in [IB]. However, the model of [IB] is based on analogy with the Vicsek IBM, 
not on an actual derivation from it. In [IB], the techniques used to show the emergence of 
a phase transition are rather complex: they are based on stability analysis in the linear 
case and renormalization group theory in the nonlinear case. In the present work, the 
model we investigate is derived from a simple IBM of collective motion combining a noisy 
Cucker-Smale consensus force and self-propulsion. The phase transition appears simply 
when analyzing the behavior of the model in the limit of a large self-propulsion force. It 
manifests itself as a change in the type of limit model at the crossing of a critical noise 
intensity. Above the critical noise, the limit model is of diffusion type, while below the 
critical noise, it is of hyperbolic type. To the best of our knowledge, the present work 
is the first instance where phase transitions in particle swarms have been evidenced in 
this way from a hydrodynamic model. A similar approach, but at the level of the kinetic 
model, can be found in [2H HH] • 

The second goal of this work is to discuss the relative merits of the Cucker-Smale and 
the Vicsek models for the derivation of hydrodynamic models of particle swarms. As men- 
tioned above, the Vicsek kinetic model imposes that kinetic velocities be of constant norm, 
while no such constraint exists in the Cucker-Smale model. Instead, a self-propulsion force 
is imposed to force the particle speed to stay close to a 'comfort' velocity. In [TO], it is 
shown that the Cucker-Smale model relaxes to the Vicsek model when the intensity of 
the self-propulsion force tends to infinity. The derivation of hydrodynamic models from 
the Vicsek kinetic model is considerably complexified by the velocity norm constraint. 
Indeed, momentum conservation is lost and the use of conserved quantities (or collision 
invariants), which is the cornerstone of the derivation of hydrodynamic models, cannot 
be implemented. In [2H], this problem has been overcome by the introduction of a new 
concept of 'Generalized Collision Invariant'. But if the hydrodynamic limit could be per- 
formed equivalently on the Cucker-Smale model, these unpleasant technicalities would be 
proven unnecessary. 

Unfortunately, performing the hydrodynamic limit on the Cucker-Smale model and 
then letting the self-propulsion force tend to infinity is not equivalent to performing the 
hydrodynamic limit on the Vicsek model. This is the second main result of the present 
work. Indeed, the type of the resulting model is the same, but the coefficients of the model 
are not the same. Specifically, the limit model has the form of a system of isothermal 
compressible Euler equations for the swarm density p and the mean velocity direction w 
(also referred to as the polarization field, see e.g. [13] )• The velocity direction oj is a vector 
of norm one. To maintain this constraint, the model includes some non-conservative terms. 
Additionally, the convective derivatives involved in the mass and momentum transport 
are not the same, a signature of a loss of Galilean invariance (see e.g. [52]). This 
model has been referred to in [25] as the 'Self-Organized Hydrodynamic (SOH)' model. 
In the present work, we show that the SOH model derived from the Cucker-Smale kinetic 
equation necessarily involves the same convective derivatives in the density and in the 
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velocity equations. Therefore, with the Cucker-Smale model, we cannot access the whole 
range of possible hydrodynamic limits that we can access with the Vicsek model. With 
the Cucker-Smale model, we only get a sub-class of these models, which limits its practical 
applications: having different convective derivatives for p and cu increases the likelihood of 
correctly reproducing emergent phenomena in swarms, such as cluster formation, waves, 
etc. 

This restriction, which is a significant disadvantage of the Cucker-Smale approach, 
can be weakened, at least partially. Indeed, in the last part of the present work, we 
include small diffusive corrections to the hydrodynamic limit of the Cucker-Smale model, 
by means of a Chapman- Enskog method (see e.g. [23] for a review). If the self-propulsion 
is taken to infinity in the resulting compressible Navier-Stokes system, then a more general 
SOH model, in particular with possibly different convective derivatives for p and u, can 
be derived. This approach is limited by the necessity to keep the diffusive corrections 
small; this limits the range of the coefficients of the SOH model which can be obtained. 
This paper shows that most SOH models can be realized in a fairly general context as 
hydrodynamic limits of either Cucker-Smale or Vicsek kinetic models. For these reasons, 
we conclude that the two approaches are somehow equivalent in the amount of technical 
work: while the Vicsek model made use of generalized collision invariants, starting from 
the Cucker-Smale model necessitates dealing with the complex diffusion terms. 

The paper is organized as follows. The problem is set up in Section [2] Then, some 
functional properties of the operators are recalled in Section [3J The hydrodynamic and 
large self-propulsion limits are derived in Section HI The diffusive corrections are dealt 
with in Section [5J Finally, a conclusion and some perspectives are drawn in Section [HJ 
Three appendices collect some of the more technical proofs. 

2 Setting of the problem 

2.1 Velocity consensus in self-propelled agent systems 

We consider a system of agents with positions Xi(t) G M. d and velocity Vi(t) G M. d , where d 
is the system dimension (in practice equal to 2 or 3), t > is the time and % G {1, . . . , N} 
is the agent's label. These agents are subject to a self-propulsion force which tends to 
restore a comfort velocity a > 0, and to a social force which drives them to the average 
velocity of the neighboring agents. Addtionally, they are subject to random velocity 
fluctuations which account for potential misperceptions and their propensity to leave the 
swarm and explore a new environment. The equations of motion are given by 




(2.1) 



with 



Ti = a 1 (v i - v^, 



J2f=i K (\ x j- x i\) 



(2.2) 
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The social force Ti is written in (12. 2p as a relaxation force towards the average velocity 
V{ in the neighborhood of particle i. The relaxation rate is a -1 (in other words, a is the 
typical time needed for agent % to align with the velocity of his neighbors). The kernel 
K, supposed spherically symmetric for the sake of simplicity, describes how the various 
partner velocities Vj are combined according to the distance of j to i. For instance, if K 
is the indicator function of the ball of radius R, it means that the agents adopt the mean 
velocity of the other agents within a distance R. The second term in the expression of 
dvi in ( 12. ip is the self-propulsion force. It takes the form of a relaxation term driving \vi\ 
towards a at rate r -1 . In other words, it takes a time r for the velocity \vi\ to relax to 
the comfort speed a. Finally, the last term is the velocity fluctuation term, where B\ are 
independent normalized Brownian processes and D > is the diffusion coefficient. The 
force Ti has been previously proposed by Cucker and Smale [22] as a model for consensus 
formation in particle swarms. A noisy version of the Cucker-Smale model is proposed in 

In the large particle limit N — > oo, by adapting the arguments in [TT], the empirical 
measure of the system 
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i=i 

where £(a* (*))(£> v ) is the Dirac delta at (xi(t),Vi(t)), can be approximated by a con- 
tinuous distribution function f(x,v,t). It solves the following Fokker-Planck equation: 



d t f + V x ■ (vf) = -V v 

with 



DAJ, (2.3) 



-r \ - J K{\x-y\)u}f(y,u,t) dydu 

r ' = ° (V '- V) - V,= SK { \ X - V \)f (v ^t)dy^- {2A) 

The left-hand side expresses particle displacement at velocity v. The right-hand side 
consists of three terms. The first one is the consensus force. The second term is the 
self-propulsion force. The last term takes into account the random velocity fluctuations. 



2.2 Scaling 

In order to understand the roles of the various terms, it is useful to introduce dimensionless 
quantities. We set x and t to be space and time units and deduce units of velocity 
= x o /to an d force Tq = ^o/^o- We assume that the range of the interaction kernel K is 
R, meaning that we can write = K{\x\/R) with K having second moment of order 

1 (i.e. J i^(|a;|)|5;| 2 dx = 0(1); we assume that K is normalized to 1, i.e. J K(\x\) dx = 1). 
We now introduce dimensionless variables x = x/xq, t = t/to, v = v/vq and the change 
of variables f(x,v,t) = Xq Vq f(xox, v v, t Q i). Finally, we introduce the dimensionless 
parameters: 

R = — , cr = — , a = — , r = — , D — —^D. (2.5) 

x t Vq t Vq 
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In this new system of coordinates, the system is written: 



d t f + V x • (vf) 



J= s = a \v f -v), 



DA v f, 



fK(^)f(y,u,t)dydw 



(2.6) 
(2.7) 



where we have dropped the tildes for the sake of clarity. Now, by fixing the relations 
between the five dimensionless parameters ( 12. 5p . we define the regime we are interested 
in. We suppose that the diffusion and social forces are simultaneously large, while the 
range of the social force tends to zero. The parameters of the self-propulsion are kept 
of order 1. More specifically, we let e <C 1 be a small parameter and we assume that 
D = 0(l/e) (large diffusion), a" 1 = 0{l/e) (large social force), R = 0(e) (small range 
of social interaction), while f" 1 = 0(1) and a = 0(1) (parameters of the social force are 
order unity). In order to highlight these scaling assumptions, we define constants D\ a% 
Rr, which are all 0(1) and such that 



D 



■D* 



a 



R = eRK 



{21 



Then, with these new notations, and dropping all 'hats' and 'sharps', we get the following 
scaled system: 



d t f e + V x -(vf e )+T- 1 V v - ((1 



-)vf 



-V v -(^f)+DA v f, (2.9) 



a-\v £ f -v) 



jK(^)uf(y,u,t)dydu} 
j x(^-)f(y,u;,t) dydu 



(2.10) 



We will investigate the limit as e — > of this system, while all other parameters (i.e. r, 
a, -D, a and R) are kept fixed. Hence, we highlight the dependence of / upon e. 

We can simplify the problem by using Taylor's expansion. At leading order, we find: 



Then, we have 
which leads to: 



Ui 



0{e% 



a 



J feu dcu 
v) + 0(e 2 ), 



(2.11) 



d t r + V x -{vf)+T- l V v - ((1 



-V, • {o-\u f , - v)f) + DA v f £ + 0(e 2 ). 



(2.12) 



6 



We will drop the 0(e 2 ) terms, since, as a first step, we consider only the leading and first 
order terms. Then, problem (I2.12p can be written as 



e d t f £ + V x ■ (vf) + r-V, • 



(2.13) 



with the collision operator Q(f): 



Q(f) 



V v - [<j-\v-u f )f + DV v f\ . 



(2.14) 



Some remarks concerning scaling ( 12. 8ft can be made. The diffusion and social forces 
are supposed of the same order of magnitude and much larger than all other forces. 
They counterbalance each other. Indeed, the social force makes the agents adopt the 
same velocity while diffusion tends to spread the velocities out. This balance results in a 
Maxwellian velocity profile (i.e. Gaussian in velocity space), as shown later on. The choice 
which is made here is to assume that the self-propulsion force is weaker. Another choice 
would have been to make the self-propulsion force as large as the social force and the 
diffusion. In this case, the balance would involve three different effects and would result 
in more complicated equilibria. This investigation is in progress [3] . The interaction range 
is supposed to tend to zero like the inverse of the interaction rate. It is no surprise that, 
at leading order, only spatially local interaction terms remain (which can be seen in ( 12. 14ft 
by the replacement of the non-local average velocity vj by the local mean velocity Uf). 
Again, other choices can be made. In [25], in the case of the Vicsek model (which is the 
limit of (I2.13P when r — > 0), it is shown that the different choice R = 0(y/e) leads to a 
different macroscopic limit when e — > 0. This choice takes better care of the non-local 
character of the interaction. It will be investigated in future work. 

Our plan is now to investigate the hydrodynamic limit e — > in this model. To this 
end, we first examine the properties of the collision operator Q. 

3 Properties of Q 

When e — > in (12.131) . f e formally converges to an element of the null-space of Q, i.e. 
a function / such that Q(f) = 0. In this limit, the dynamics are characterized by the 
projection of the left-hand side of (12.131) onto the space orthogonal to the range of Q. This 
space is spanned by the so-called 'collision invariants'. In this section, we successively 
determine the null-space of Q and its collision invariants. 



We first define the Maxwellian with mean velocity u G M. d and temperature T = oD > 



3.1 Null-Space 



as follows: 




(3.15) 
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Note that M u satisfies J M u dv = 1 and J M u v dv = u. 

To proceed, we need to determine a functional setting. Let u £ M. d and define a 
weighted L 2 -space H u such that 

and a weighted if 1 -space V u such that 

K := {<P e H u : / |V„0| 2 M u cfo < +cx)} 



with the associated norms: 



l M u dv, 



Ik 



|V„0| 2 M u ^, 



2 ^ = I^IL + 



(3.16) 



and inner products (•,•)#■„, (-, -)v u , and ((-,-)) 14 respectively. 

Lemma 3.1. (i)The operator Q given by f)2.14p can be reformulated as: 



Q(f) = DV V 



M u yJ—) 

1 v M ' 
iVlu i 



(3.17) 



(ii) For any function f{v) such that f /M u G V Uf and for any function g £ V Uf , we have: 

j Q(f)gdv = -D J M Uf V v (j^)-V v gdv. (3.18) 
In particular, we have: 



/W)(£) *v = -dJm u , 



dv. 



(3.19) 



(Hi) The null-space KerQ = {f(v) \ f/M u . £ V Uf and Q(f) = 0} is given by: 

KerQ = {pM u \ p > 0, u £ R d }. (3.20) 

The proof of this lemma is postponed to Appendix El An element pM u of Ker(Q) is called 
a local thermodynamic equilibrium with density p and mean velocity u. 

3.2 Collision invariants 

Definition 3.2. A function ip(v) is said to be a collision invariant (CI) if and only if 

Q{f)1>dv = 0, 



for every f such that f/M Uf £ V u , and ip £ V Uf . The set of CI's is denoted by C. It is a 
vector space. 
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We have the: 

Proposition 3.3. We have C = Span{l, v%, . . . , v^}. In other words, vj is a CI if and 
only if there exists a £l and b G M. d such that ip(v) = a + b ■ v . 

The proof of this proposition is again postponed to Appendix lAl 



4 Hydrodynamic limit and fast relaxation 

The goal of this section is to investigate the formal limit e — > in ( 12.131) and to examine 
some of the properties of the limit system relative to the propulsion force. More precisely, 
we exhibit a phase transition when the intensity of the velocity fluctuations crosses a 
certain threshold dependent on the magnitude of the propulsion velocity. 



4.1 Hydrodynamic limit 

The goal of this section is to prove the following formal theorem: 

Theorem 4.1. Let f e be the solution of equation ^2.13\) associated to an initial datum 
fi(x, v). We suppose that fi is independent of e for simplicity. We assume that solutions 
of Ii2.13\) exist on any time interval [0, 7~]. Assume that f £ — > f° as e — > as smoothly 
as needed, which means in particular that derivatives of f e converge to the corresponding 
derivatives of f°. Then, there exist two functions p(x,t) > and u(x,t) 6 M d such that 

f°(x,v,t) = p(x,t)M u(x , t) (v), V(x,v,t) G R 2d x [0,71. (4-1) 

Furthermore, p and u satisfy the following system of isothermal compressible Euler equa- 
tions with relaxation: 

8 t p + V x . • (pu) = 0, (4.2) 
d t (pu) + V x • (pu <g> u) + TV xP = -lp M ( l '" |2 + ( ^ + 2)T - 1) , (4.3) 
associated with initial data (pi, ui) such that 

Pi = fi dv, piuj = fiv dv. 



Proof. First, we note from (1243]) that Q(f £ ) = 0(e). Therefore Q(f°) = 0, which, 
because of (I3.20p . implies that f° is of the form (14.11) . 

Next, using the CI's given by Proposition 13.31 we multiply ( 12 . 1 3[) successively by 1 
and v and use the fact that the right-hand side vanishes upon integration. Then, we get 
the following conservation relations, which are valid for any s: 

d t p £ + V x -f = 0, (4.4) 
d t f + V, • E £ = r- Y, (4.5) 
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with p £ , j £ , £ £ , the density, flux and pressure tensor associated to / £ , given by: 

p £ = f f £ dv, f = [ fv dv, £ £ = [ f (v <g> v) dv, (4.6) 



and right-hand side q e given by 



\v 



2 



q e = / fvil-^dv. (4.7) 



a 2 



Now, letting e — > 0, we can express j = lim j £ , E = lim£ £ and q = limg £ as functions 
of p and u: 



j = pM u v dv = pu, 



E = / pM u (u ® v) dv = p(u ®u) + pTId, 



\v 



2 



u\ 2 + (d + 2)T 



q= pM u v(l -)dv = pu( 1). 

er a z 

Inserting these expressions into the conservation equations leads to (14. 2p . (I4.3p . The 
statement about the initial conditions is obvious. ■ 

Remark 4.1. Dividing (14. 3 p by p and using ( 14.2ft . the momentum conservation equation 
can be written in n on- conservative form: 

(d t + u- V x > + T— ^ = —^u(\u\ 2 - (a 2 -(d+ 2)T)) . (4.8) 



This model is nothing but the Euler system of isothermal compressible gas dynamics with 
a forcing term. The mass conservation equation (14 ,2p (also known as the continuity equa- 
tion) has a standard form. The momentum balance equation (14. 3 \ involves an isothermal 
pressure Tp on the left-hand side and a self-propulsion force on the right-hand side. The 
temperature T = aD is proportional to the ratio of the intensities of the velocity fluc- 
tuations D and of the social force cr _1 . The self-propulsion force takes the form of a 
relaxation of the fluid velocity to a comfort fluid velocity ^ a? — (d + 2)T. In comparison 
to the force acting on individual particles, the force acting on the fluid involves a term 
depending on the temperature. The temperature being a truly macroscopic quantity, 
this term cannot have any counterpart at the particle level. Here, at the fluid level, the 
comfort fluid velocity becomes a pure imaginary number when T is larger than a critical 
temperature T c = a 2 jid + 2). When the temperature crosses T c , a phase transition oc- 
curs. This phase transition is studied in the next section when the intensity r _1 of the 
self-propulsion force is taken to infinity, a limit which we refer to as the 'fast relaxation 
limit'. 
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4.2 Fast relaxation limit in the hydrodynamic model 

We recall that, according to our notation, r is the ratio of the physical relaxation time 
to the time unit to (see Section 12. 2\\ and is a dimensionless parameter. The goal of this 
section is to investigate the limit r — > in the hydrodynamic model (14.21) . (14. 3p . We 
denote by (p r , u T ) the solution for finite r and assume that its limit (p, u) as r — > exists 
and is as smooth as needed. We will use the non-conservative form of the model, which 
we recall here for the sake of convenience: 

d tP T + V, • (p T u T ) = 0, (4.9) 
(d t + u T • V> T + TV X hip- = —u T C ' l 2 " !)• ( 4 - 10 ) 

Letting r — > formally in (!4.1Up leads to |u| 2 + (d + 2)T — a 2 = 0. Therefore, there are 
two cases according to whether the quantity (d + 2)T — a 2 is positive or negative, i.e. 
according to the position of T with respect to the critical temperature T c defined by 

a 2 

T c = -i—. (4.11) 
c d+2 v ; 

4.2.1 Case T > T c (large noise) 

We let (d + 2)T - a 2 = (d + 2)(T - T c ) := s 2 > 0. The constant s only depends on the 
problem data and not on the solution. In this case, equation (I4.10p becomes 

(d t + u T - V x )u T + TVJnp T = —L u T (\u T \ 2 + s 2 ). (4.12) 

TCL Z 

To examine the limit r — > 0, we need to find the equilibria of the right-hand side of (14. 12ft . 
i.e. the solutions u of u(\u\ 2 + s 2 ) = 0. Obviously, the only solution is u = 0. Additionally, 
at least in the spatially homogeneous setting, this is a stable solution. Indeed, the unique 
solution of 

— = -u{\u\ +s), w(0)=w , 

satisfies u(t) — > as t — > oo. So, in the spatially non-homogeneous case, we formally have 
u T — > as r — > 0. Therefore, the formal limit of (14. 2p . f)4.3p gives 

d t p = 0, n = 0. 

In order to get a more precise description of the limit r — > 0, we need to rescale time and 
velocity. 

With this aim, we let t' = rt and u T (x,t) = ru T (x,t'), p T (x,t) = p T (x,t'). Inserting 
this into (14.21) . (14.31) . we find (dropping the tildes): 

d t p T + V x ■ (p T u T ) = 0, (4.13) 

r 2 (d t + u T ■ V x )u T + TV X lnp r = -\u t {t 2 \u t \ 2 + s 2 ). (4.14) 

The behavior of this system when r — > is that of a diffusion. More precisely, we state: 
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Proposition 4.2. Assume that the solution (p T ,u T ) of system M-lcfy , Ij4-14\ ) is smooth 
and converges smoothly towards a pair (p,u) as r — >• 0. Then, p satisfies the following 
diffusion equation: 

TT 

d t p-D dlff A xP = 0, Z^ iif =__i_. (4.15) 

and u = -D di ffV x \np. 

Proof. Taking r -»• in equation (14.1411 . we get u = -^fV x lnp = -jr^rV x \np. The 
limit model therefore follows from the continuity equation ( 14 . 1 3 1) and is given by (14. 15 j) . ■ 

Remark 4.2. We note that D^g — > oo when T — > T c . Therefore, the time variation of 
p is faster as T — T c gets smaller. We can rescale time to a faster time-scale by setting 
t' = D ^ff t = jP-i t. Then, in the rescaled variables, the diffusion equation ([^. i<5| ) becomes 
independent of T — T c and gives the standard heat equation d t p — A x p = ( omitting 
the primes for simplicity) . 

4.2.2 Case T <T C (small noise) 

Let us now consider the dynamics of the system in the case of smaller noise. In this 
section, we aim to prove the following: 

Proposition 4.3. LetT < T c and define c 2 = a 2 -(d+2)T = (d+2)(T c -T) > 0. Assume 
that the solution (p T ,u T ) of system d^.i0| ) is smooth and converges smoothly towards 

a pair (p,u) as r — > 0. Assume additionally that (p T ,u T ) is not identically equal to {po,0) 
where po constant in both space and time. Then, u = cu, where u G S^" 1 and the pair 
(p, oj) satisfies the following system: 

d t p + V x ■ (cpu) = 0, (4.16) 

(d t + cco-V x )co + -P(V x lnp) = 0, (4.17) 

c 

where P = Id — u> (g) oj is the orthogonal projection onto the hyperplane orthogonal to u>. 
Remark 4.3. We easily see, taking the scalar product of j4-ll ) with u, that 

{d t + cuj-V x )\u\ 2 = 0. 
This implies that \ui{x,t) \ = 1 for all time provided that \u(x,0)\ = 1 initially. 

Proof. Equation (14. 1 Oj) is now written as: 

(d t + u T ■ V> r + TV X lnp r = \u T (c 2 - \u T \ 2 ) . (4.18) 

TGT 

We now look for the equilibria of the right-hand side of ( 14 . 1 8 [1 . namely the solutions of 

u ( c 2 - |m| 2 ) = 0. 
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There are two sets of equilibria. The first set reduces to the single point u — 0. The second 
set is the sphere \u\ 2 = c 2 . We now show that in the spatially homogeneous setting, the 
first equilibrium is unstable, while the second class of equilibria is orbitally stable. Indeed, 
let us consider the solution u(t) of the differential equation 

diUj / o i i o \ / \ 

— = u(c - \u\ ), u(0) = Wo- 
ItS solution can be analytically given by 



u{t)\ 2 = 





\u 


2 c 2 


(c 2 - 




2 )e 


,-2cH + 


U 


2 



and is such that \u(t)\ c for all initial data Uq except Uq = 0. This shows that uq = 
is an unstable equilibrium. On the other hand, if one perturbs an equilibrium \uq\ 2 = c 2 
by a small amount, the solution will relax to an element of the circle \u\ 2 = c 2 (where u 
may be different from uq). 

Now, we let r — > in the non spatially homogeneous system (14.161) . (14.181) . Unless u T 
is identically zero, which can only occur for a uniform density po, u T converges (at least 
formally) towards one of the stable equilibria. Therefore, we have u T — > u = cu with 
\u\ = 1. In order to find the equation satisfied by u, we introduce a polar decomposition 
of the solution u T . We write u T = c T u T , with c r = \u T \ and u T = u T /c T . Let P T the 
orthogonal projection of M. d onto the hyperplane orthogonal to u T . The projection P T can 
be written tensorwise as P T = Id — u T <E> oo T . Inserting the polar decomposition of u T in 
pg| leads to: 

co T (d t + u T ■ V x )c T + c T {d t + u T ■ V x )u T + TV X lnp r = -^u T (c 2 - (c r ) 2 ) . (4.19) 

We note that, because \u T \ = 1 and the operator d t + u T ■ V x is a derivative, the vector 
(d t + u T ■ V x )u T is orthogonal to u T . Furthermore, the first term of the left-hand side and 
the right-hand side are parallel to u T . Therefore, applying P T to the second term of the 
left-hand side leaves it unchanged, while it cancels the first term of the left-hand side and 
the right-hand side. Consequently, applying P T to (14.191) leads to: 

c T (d t + u T -V x )uj t + TP T (V X In p T ) =0. (4.20) 

Now, taking the limit r — > formally leads to (14.171) . ■ 

System (I4.16p . (14. 17ft has been referred to in the literature as the Self-Organized Hydrody- 
namic (SOH) system [22J[25]. It has the form of a compressible gas dynamics system with 
isothermal equation-of-state and geometric constraint \u\ = 1. The projection operator 
P which multiplies the pressure term V x In p maintains the constraint over the course of 
time (see Remark 14. 3ft . It results in a non- conservative term (since P depends on u). 
This type of system has been derived for the first time in [26] and has been shown to be 
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hyperbolic. Beyond this result, the mathematical study of such systems is in its infancy. 
A local existence result is given in [23] and some special solutions are given in [H]. 

When T -4 T c , we have c — > and therefore, p becomes constant. To find a non- 
trivial dynamics, we must rescale time to a slower time-scale. By the time rescaling 
t' = ct, system (I4.16p . (14.171) can be written (dropping the primes): 

dtp + V x ■ (pu) = 0, 

(9 ' + W ' V + (d + 2) T (Tc-T) P ^' = °< 

The parameter (T c — T)/T plays the role of the squared Mach-number in the standard 
compressible Euler system. Therefore, the limit T — > T c , is similar to a small Mach- 
number limit. When T A T c , we formally get p — > po where po is a constant in space. 
With appropriate boundary conditions we can assume that po is also independent of time. 
Then, the limit system as T 4 T c is written as follows: 

V x ■ u) = 0, 

(dt + cu- V x )cu + PV,tt = 0, 
where 7r is a hydrodystatic pressure defined by 

This system has been already proposed in [27] for modeling gregariousness. It is a system 
of incompressible Euler equations subject to the geometric constraint \oj\ = 1. 

The behavior of the (appropriately rescaled) SOH system (14. 16ft . ( 14. 17ft when t4t c 
is very different from the behavior of the (appropriately rescaled) diffusion equation (14 . 1 5 f) 
when T A T c . The former is given by the incompressible Euler equations with geometric 
constraint \tu\ — 1 while the latter is the standard diffusion equation (see Remark 14.21 ). 
This drastic change of type is the signature of a phase transition further elaborated on in 
Section ELl 



4.2.3 Comments 

The fast relaxation limit of the compressible Euler system with self-propulsion (14. 2 p , (14. 3 p 
exhibits two different regimes for temperatures below or above the critical temperature 
T c . Above T c , the behavior of the system is that of a diffusion, while below T c , the system 
obeys a hyperbolic system, the Self-Organized Hydrodynamic (SOH) model (I4.16p . (I4.17p . 
Even when the temperature is close to the critical temperature, it is not possible to match 
the two types of models in a smooth way, as noticed in the previous section. This abrupt 
change in the type of the model as the temperature crosses a threshold is a manifestation 
of a phase transition. Since T c depends on a, the phase transition originates from the 
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Kinetic model (1235]) 



t-s-o no] 



Kinetic model 
with norm constraint on v 



e-s-0 (section [47TJ 



e->0 [26] 



Compressible Euler 
with self propulsion 



r^o (section l4T2)i Self- Organized Hydrodynamics 



HED, (SUD 



Figure 1: Two strategies. First strategy: first take the relaxation limit r — > (top hori- 
zontal arrow) and then pass to the hydrodynamic limit e — > in the resulting constrained 
kinetic model (right vertical arrow). This is done in [26] and [10]. Second strategy: first 
pass to the hydrodynamic limit e — > in the original kinetic model (left vertical arrow) 
and then pass to the fast relaxation limit r — > in the resulting Euler system (bottom 
horizontal arrow). This is what is done here. 



self-propulsion force. Indeed, when a = (i.e. no self-propulsion), T c = and there is no 
phase transition: the system is in the diffusion regime in all circumstances. 

Phase transitions in self-propelled particle systems have already been evidenced [S3] 
and an abundant literature has been devoted to them (see the review [M])- Phase transi- 
tion in hydrodynamic models of self-propelled particles have first been studied in [IB] (see 
also the review [19]). However, the equations proposed in [18] are more complicated than 
the ones seen here. They are derived solely on heuristic principles and invariance con- 
siderations. Their analysis is based on a combination of linear and nonlinear techniques. 
There is no link to the underlying particle models. The link between hydrodynamic and 
particle models of self-propeled particle systems has been made in jTJ IB], but for binary 
interaction mechanisms instead of the mean-field interaction considered here. Here, the 
hydrodynamic model is derived from the underlying particle dynamics and is much sim- 
pler: it merely consists of the isothermal compressible Euler model complemented with 
the self-propulsion force. The phase transition manifests itself in the change of type of 
the PDE which describes the system under large self-propulsion. 

The SOH model (I4.16p . (I4.17P has previously been derived in [26] from a system of 
self-propelled particles which have constant and uniform velocity. In [ID] , it has been 
shown that the kinetic model with the velocity norm constraint of [2S] is the the fast 
relaxation limit r — > of (I2.13p . A natural question is then whether the imposition of the 
norm constraint at the particle and kinetic levels as in [26] is necessary or useful. Indeed, 
there are now two ways of deriving the SOH model from ( I2.13p . which are summarized in 
Fig. [TJ The first way is to follow the top horizontal and right vertical arrows successively. 
This is what is done in [26] and [10]. The second way is to follows the left vertical and 
bottom horizontal arrows successively. This is what is done here. 

There is however one noticeable difference between the two strategies. In [5E] , the 
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system is written: 



dtp + Vs • (cxpu) = 0, 

(d t + c 2 u ■ V x )u + 8P(V X hip) = 0, 



(4.21) 
(4.22) 




where the constants c\ and c 2 are such that < c 2 < c\. In [33], it has even been shown 

that taking into account an anisotropic vision, (ci,c 2 ) can be made arbitrary. We would 
recover (I4.16p . (14.171) if we could set c 2 = C\ = c and 5 = T/c. However, it is not possible 
since, in [26J, c 2 ^ c\. Therefore, the diagram in Fig. [T]is not commutative. This shows 
how important it is to enforce the norm constraint at the kinetic level, in spite of the 
induced difficulties. Indeed, such constrained kinetic models do not exhibit momentum 
conservation and lack classical collision invariants such as those derived in Section 13.21 
This makes the derivation of hydrodynamic models considerably more complex. In , 
new concepts have been developed to bypass these difficulties. The fact that c\ ^ c 2 is a 
direct consequence of these features. 

The question whether c\ = c 2 or not has important consequences. Indeed, the SOH 
system is Galilean invariant if and only if c\ = c 2 (which is the case here). However, in 
the generic case C\ ^ c 2 , the SOH system is not Galilean invariant. This fact reflects a 
key feature of collective motion: the anisotropy of information flow. This is evidenced 
in car traffic, where perturbations of velocities (typically moving with speed c 2 < c\) 
propagate upstream the flow (whose speed is ci). This phenomenon is a consequence 
of the fact that information propagates upstream from drivers ahead to drivers behind. 
When Ci = c 2 , this property is lost and information spreads in an isotropic way just 
like usual gas dynamics. In such a case, the model is unable to correctly reproduce the 
complex emerging patterns (such as congestions and waves in car traffic). We refer to [52] 
for a discussion of the loss of Galilean invariance in biological swarms. 

The approach of [26] does not lead to phase transitions: whatever the values of the 
temperature T or self-propulsion velocity a, the limit system is of hydrodynamic type. 
The very large intensity of the self propulsion force which is needed to perform the r — > 
limit first (top horizontal arrow in diagram [T]) prevents any diffusion regime to establish. 
However, the emergence of a phase transition is a well-established experimental fact (see 
review in [M])- The inability of [26] to produce phase transitions can be seen as a major 
drawback. Fortunately, it has been proved in [21] that adding a dependence of the social 
force intensity a" 1 upon the flux j restores the phase transition. It also brings density- 
dependent phase transition. More precisely, the parameter controling the phase transition 
is the ratio T/p. A natural question is to explore similar features here. This question will 
be investigated in future work. 

Here, we choose a different direction: we explore whether the inclusion of diffusion 
terms can cure the deficiency of the SOH models derived here (i.e. the fact that C\ = c 2 ). 
This is the goal of the next section. 
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5 Diffusive corrections to the hydrodynamic model 



5.1 Setting of the problem 

In this section, we derive 0(e) diffusive corrections to the compressible Euler model with 
self propulsion (14. 2p . (14. 3p . We show that these diffusive corrections lead to a system of 
isothermal compressible Navier-Stokes equations including self-propulsion. In a second 
step, we perform the fast relaxation limit r — > in the resulting Navier-Stokes model and 
compare it to the SOH model obtained in Proposition 14.31 
We now set up the problem. We define 



Pf = Jfdv, j f = J fvdv = p f Uf. (5.23) 

We assume that the observation kernel if (|£|) is such that J K(\£\) d£ — 1 and we denote 
by k > the second moment of K, i.e. 

\ f K{\t\){i®i)dii = kU. 

For instance, if K is the indicator function of the ball of radius 1, k — 2 lf(rf+2) ' wnere 
\§> d ~ 1 \ is the (d — l)-dimensional measure of E> d ~ l . In the cases of d = 2 and d — 3, we 
respectively get k = n/8 and k = 2it/VS. 

We first give the expansion of J^j up to the fourth order in e. 



Lemma 5.1. We have: 



=u f + e 2 u) + 0(e 4 ), (5.24) 
jr^ = a ~ l {u f -v)+ e 2 a~ 1 u 1 f + 0(e 4 ), (5.25) 



where 



"f 



) = {PfA x jf - j/AxPf) = k R (A x u f + 2(V x .lnp/ • V x )u f ), (5.26) 
Pf 



and k R = kR 2 



We give the proof of this very simple lemma in Appendix [Bj The 0(e 2 ) correction to the 
mean velocity Uf takes into account the non-local character of the average 1 12. 10p . This 
correction involves gradients of the local density and flux. They quantify how information 
spreads due to the fact that the agents observe their environment over a certain spatial 
extent. Because the observation is supposed isotropic, the correction only involves 0(e 2 ) 
terms and second order derivatives. In the case of non-isotropic observation kernels, 0(e) 
corrections involving first order gradients would be obtained. The study of this effect is 
postponed to future work. 
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Now, up to terms of order e 4 (which are dropped), equation (I2.9P is written: 



dtf + V x ■ (vf £ ) + r^V, ■ ((1 - ^)vf £ 

+ e 2 a~ l V v -(u) £ f £ ) 



Q(f £ 



(5.27) 



where again, Q is given by (I2.14|) . Now, we look for a fluid model approximating ( I5.27P 
which includes the 0(e) terms. We will adopt Chapman-Enskog's method consisting in 
writing the 0(e) terms in terms of spatial derivatives only (see e.g. the review in [23]). 
The model is introduced and discussed in the next section. 



5.2 Compressible Navier-Stokes equations with self-propulsion 

The compressible Navier-Stokes system with self-propulsion is established in the following: 

Theorem 5.2. Let f £ be the solution of associated to a given initial condition fj, 

which is supposed independent of e for simplicity. Let (p/e, pfeUfe) be the moments of f e 
defined by i5.2cti) . Then, we can formally write (pfe, pjeUfe) = (p £ , p £ u £ ) + 0(e 2 ), provided 
that (p £ ,p £ u £ ) satisfy the following set of compressible Navier-Stokes-like equations (where 
we drop the superscripts e upon p £ and u £ for the sake of clarity): 

dtp + V x ■ (pu) = 0, (5.28) 

1 Id 2 

d t (pu) + X £ V X ■ (pu <S>u) + V x ti £ (p, \u\) = — -pu(^r ~ X s ) 

t £ a z 

+ e(pV x • (pS(u)) H -pA x u) H (V x p ■ V x )u 

a a 

e\ \u\ 2 
+ y P [(Vx • u)u + V x (i^-) + (u ■ V x )u] , (5.29) 



where 



\ £ = l-e\, A=^ r , (5.30) 
ra 2 

7r e (p, = Tp- £vr(p, \u\), 7r(p, |«|) = ^p{(d + 2)T - a 2 + |u| 2 }, (5.31) 

x e = l-f(d + 2)-(d + 2)^(l-f(d + 4)) ^ (5 ^ 



2 



(d 



i = I(l_ e ^(d + 8 )), (5.33) 
p = aT, E(u) = ^(V x u + (V,m) t ). (5.34) 



The proof of this result is fairly technical. It is given in Appendix O The coefficient p is 
the fluid viscosity, while ii £ (p,u) is a velocity-dependent pressure. 
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We compare this Navier-Stokes model to the hydrodynamic model (j4.2p . (14. 3p and 
provide a physical interpretation of the 0(e) correction terms. The mass conservation 
equation (I5.28P is unchanged compared to (14.21) . as it should be. Indeed, the particle 
density is still a conservative variable transported by the fluid velocity. By contrast, there 
is a wealth of new terms in the momentum balance equation (I5.29p . We can identify the 
only term which come from classical fluid viscosity: this is epV x ■ (pS(u)). The present 
viscous term is different from that appearing in the full Navier-Stokes system, where 
the temperature is determined by the energy balance equation. Here, the viscous term 
involves the symmetrized velocity gradient tensor S(u) instead of the rate of strain tensor 
S(u) = S(u) — ^(V x • u)Id. This difference is solely due to the isothermal character of the 
model and not to self-propulsion. 

Then, there are two terms which come from the non-locality of the interaction and 
which are proportional to /cr. The first one contributes to adding more viscosity and is 
equal to e^pA x u. The second one contributes to convecting the velocity in the direction 
of the gradient of p. It is written ^^(V x p ■ V ' x )u. 

All other 0(e) terms originate from self-propulsion. The self-propulsion first con- 
tributes to a similar relaxation source term as in (jOjl : the term -±pu(^jf - x £ )- How- 
ever, both the relaxation rate \ and the bulk comfort velocity y/x^a towards which 
this source term is relaxing are different from those of (I4.3p . respectively equal to - and 
\J a? — (d + 2)T. We notice that the former are O(e) corrections of the latter, as they 
should be. However, there is a second source term, the term 

e\ \u\ 2 

yp[(V,-«)« + V4^) + (tt-V>], 

which describes a force comprising three terms. The first term is just a friction propor- 
tional to the compressibility V x • u. The second term is a force acting in the direction of 
the gradient of \u\ 2 . Finally, the last term is a force acting in the direction of (u ■ V x )u. 

The self-propulsion force also induces some changes in the convection terms. The first 
one is a change in the convection velocity of the momentum, which is not just u but \ £ u. 
The coefficient A e is close to one but not exactly equal to one, the difference being 0{e). 
This feature makes the model closer to the generic SOH model ( 14.211) . (14. 22ft . which has 
different convection velocities c\ and C2 for the density and velocity. The last contribution 
of the self-propulsion is a modification of the pressure. The isothermal pressure pT is 
complemented by an 0(e) pressure correction eir(p, \u\) which depends on the norm of 
the velocity. This correction is positive or negative according to whether the speed is larger 
or smaller than the bulk comfort velocity \J a? — (d + 2)T of the non-viscous model. 

We now discuss this model in view of the model proposed in [IBJ H9], which is consid- 
ered as the 'paradigmatic' model of hydrodynamic type for flocking. In this model, the 
mass conservation equation is the same as (I5.28P but the momentum balance equation 
takes the following form (in the absence of extremal force): 

d t u + (u ■ V x )u = au - (3\u\ 2 u - V X P + D L V X (V X ■ u) 

+DiA x u + D 2 (uV x -u) 2 u, (5.35) 
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where a, (3, D L , D 1 and D 2 are positive coefficients and P = P(p) is the pressure, 
which depends nonlinearly of p. This model has been derived on the basis of invariance 
considerations. In order to compare our model with (15. 35ft . we use the mass conservation 
equation ( 15. 28ft to write the momentum balance equation (15.291) as follows: 

d t u + (l - £y) (u ■ V x )u -e\(u- V x . hip) u - £y (V x -u)u-e^ V x (^j-) 

1 \u\ 2 21cr 
= -pu{-^- - X s ) ~ V a .7r e (p, \u\) + e(p-\ )(V x ln/9 ■ V x )u 

+e(p^ )A x u + ep(V x u)V x \np + epV x {V x -u). (5.36) 

There is one term in (15.351) which is missing from f l5.36[) : the anisotropic velocity diffusion 
D 2 {u- V x ) 2 u. By contrast, there are many terms appearing in (I5.36P which are not present 
in (15.351) : the third, fourth and fifth terms of the left-hand side and the third and fifth 
terms of the right-hand side. Additionally, among the terms which are common to both 
formulas (namely the first and second terms of the left-hand side and the first, second, 
fourth and sixth terms of the right-hand side of (15.361) ) , some of them assume different 
forms. Indeed, the second term (the convection term (u ■ V x )u) is multiplied by the 
constant (l — e^-) less than one in (15. 36 p . while this coefficient is exactly one in (I5.35p . 
This difference is significant, in view of the previous discussion about mass and momentum 
convection velocities in swarming systems. Another difference is in the pressure term 
(second term at the right-hand side of (I5.36P ). In our model, the pressure depends on 
both the density and the norm of the velocity, while it depens on the density only in (I5.35p . 
By contrast, the dependence upon the density is linear in our case, while it is nonlinear in 
(I5.35p . This discussion illustrates that phenomenological models can differ significantly 
from first principle models when complex phenomena such as swarming behavior are 
concerned. The question whether these differences lead to perceivable changes in the 
qualitative behavior of the solution has not yet been investigated. 

The mathematical properties of this system (such as e.g. the stability of this system 
for small perturbations of a homogeneous state) will be studied in future work. In this 
work, we investigate the fast relaxation limit r — > 0. This is the goal of the next section. 

5.3 The fast relaxation limit in the compressible Navier-Stokes 
equations with self-propulsion 

In this section, we examine the limit r — > in the Navier-Stokes system (I5.28p . (I5.29|) . 
Since this system was derived under the assumption that e is small, e needs to tend to 
at least as fast as r tends to 0. Here, we decide to make e and r proportional. This 
is the borderline case, because, to be consistent, we should have linked e and r in this 
way already at the kinetic level. However, as pointed out earlier, the analysis of this 
scaling limit is more complex and is still under scrutiny [3]. Our conjecture is that the 
kind of model we get is the same in both limits. Therefore, investigating it at the level of 
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the Navier-Stokes system is a preparation before performing the limit directly from the 
kinetic level. 

Since A is proportional to 1/r, we decide to relate e and A in such a way that eX is a 
constant a, i.e. 

a = eX. (5.37) 

We need to keep in mind that, strictly speaking, a must be 1 otherwise the derivation 
of the Navier-Stokes model in the previous section loses its valididty. We decide to express 
e as a function of r. We can write 

aa 2 

e = n a r, K a = ——. (5.38) 
Zal 

We also define 

A a ;=A^ = l-a, (5.39) 
7r a (p, := 7rf(p, \u\) = Tp- |p{(rf + 2)T - a 2 + |u| 2 }, (5.40) 

1 1 & t 1 d + 8 {ZA1\ 

— = — = —, £ a = 1 —a. (5.41) 

r Q t z 

We can introduce an a-dependent critical temperature T c (a) by: 

1 - 

T e (a)=T e (0)- jfj-, (5.42) 

1 - -±-a 

where T c (0) is the critical temperature ( 14. lip . Then, we can write the square of the 
comfort velocity in the Navier-Stokes model c\(a) as 

1 - 

cl(a) = a 2 X - =(d + 2) jL- (T e (a) - T). (5.43) 

Since a must be small, we can limit its range in such a way that £ a > and that the 
relaxation time r a remains positive. Therefore, we have a G [0, ^rg]- In this interval, 
T c (a) is an increasing function of a with values in [T c (0), |T C (0)]. We also assume that 
the temperature is below the critical temperature T < T c (a) in such a way that cf (a) > 0. 

We finally denote by (p T ,u T ) the solution of the Navier-Stokes system (15.281) . ( 15.291) . 
We re-write the system in the new notation: 

d tP T + V, ■ (p T u T ) = 0, (5.44) 
d t (p T u T ) + X a V x ■ (p T u T (g> u T ) + V x 7r a (p\ \u r \) = -- ^p T u T (\u\ 2 - c\{a)) 



t a 



+ r 



— ( (pV x ■ (p T £(u T )) + ^ P T A x u T ) + ^(V xP r ■ V x )v7 



a r ,_ ^ ~ _ ,\u T[ 



2 p[(V, • u T )u T + V.(^) + (u T ■ V x )u T ], (5.45) 

Now, we formally let r — > in this system, keeping all other parameters fixed, and in 
particular a. We get the following: 
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Proposition 5.3. Let a G [0, ^Tg] and let T < T c (a). Assume that the solution (p T ,u T ) 
of system ft5.44\ )> ( |5.^<5| j is smooth and converges smoothly towards a pair (p, u) as r — > 0. 
Assume additionally that (p T ,u T ) is not identically equal to (po> 0) where po is constant in 
both space and time. Then, u = ci(a)u, where u G S d_1 and the pair (p,to) satisfies the 
following system: 

d t p + V x -(c 1 {a)pu)=0, (5.46) 
(dt + c 2 (a)u ■ V x )u + 6 a P{V x \np) = 0, (5.47) 

where P = Id — u (g) uj is the orthogonal projection onto the hyperplane orthogonal to u 
and 

3 

Pa (a) = (l - -a)ci(a), (5.48) 

*- = ^)> r « = T^{( 1 + ^) r -H' (549) 

Remark 5.1. For small a, up to terms of order a 2 , we have: 

cx{a) = ci(0) + a(d + 2)(3T C (0) - 2T) > Cl (0), 
c 2 {a) = ci(0) + a^^(3T c (0) - T) > c 2 (0). 

Therefore, the convection speeds are larger when first order corrections are included. Ad- 
ditionally, it is easy to prove that c\{a) is an increasing function of a on its interval of 
definition. 

Proof. The proof is similar to that of Proposition 15.31 and is only sketched. First, 
observe that the second line of the right-hand side of (I5.45p . which is proportional to r, 
simply vanishes in the limit. Since all the diffusion terms in the Navier-Stokes system are 
contained on this line, there is no diffusion in the limit system. Second, since \u T \ — > cx(a), 
which is a constant, the second term of the third line of (15.451) also vanishes in the limit. 
For the same reason, the pressure n a (p T , \u T \) —> vr Q (p, Ci(a)) = T a p. Therefore, we recover 
an isothermal pressure equation-of-state in the limit. Then, since (15.491) is obtained by 
projecting (15.451) onto the hyperplane normal to u T , the first term of the third line of 
(I5.45p . which is parallel to u T , also vanishes in the limit. Finally, the last term of the 
third line of (I5.45P combines with the second term at the left-hand side of (I5.45P and 
yields the second term of the left-hand side of (15.471) . Indeed, we get 

P T [A a V x • {p T u T ® u T ) - -p{u r ■ V x )u T ] -> (A Q - -)pc\{a){uj ■ V x )co, 

where P T is defined like in the proof of Proposition 15.31 Therefore, c 2 (a) = (A Q — |)ci(a) 
is given by (I5.48H . The remaining details are left to the reader. ■ 
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Now, system flOgj) , flo~4Tj) is in the form of the generic SOH model (H~2T]) . f l4~2"2"]) with 
c 2 < ci, like in We conclude that the inclusion of diffusive terms in the hydro dynamic 
model before passing to the fast relaxation limit was a successful approach. It has cured 
the deficiency of the SOH model derived in Section 14.2.21 In some sense, the purely 
hydrodynamic model obtained at Theorem 14. II is too simple to describe the full complexity 
of the system but the inclusion of diffusive terms is enough to restore the adequate level 
of complexity. 

It is intriguing that the only diffusive corrections that are kept in the fast relaxation 
limit are those coming from the self-propulsion force. It would be interesting to compute 
the 0(t) corrections to this model. Then, the other diffusive terms (those arising from 
viscosity and non-locality of the interaction) would appear. The resulting model would 
then have more relevance as an approximation of the original kinetic model. It should 
also be compared to the diffusive SOH model obtained in [28], which is quite complex. It 
would be instructive to see if the present approach could help get a cleaner model. 

6 Conclusion and perspectives 

In this work, we have provided evidence of a phase transition from disordered to ordered 
motion in a hydrodynamic model of socially interacting agents with self-propulsion. The 
model we have investigated has been derived from a particle system combining a noisy 
Cucker-Smale consensus force and self-propulsion. We have shown that the phase transi- 
tion appears in the limit of a large self-propulsion force and manifests itself as a change 
of type of the limit model (from hyperbolic to diffusive) at the crossing of a critical noise 
intensity. We have also shown that, in the hyperbolic regime, the resulting SOH (self- 
Organized Hydrodynamics) model suffers from unnecessary restrictions on the range of 
its coefficients. To remove these restrictions, we have computed diffusive correction to the 
model. With these diffusive corrections, the restrictions on the SOH model obtained in 
the limit of a large propulsion force disappear. 

As pointed out in the core of the work, many points deserve further elucidation. A 
first one, currently under scrutiny, consists in performing the combined hydrodynamic and 
large self-propulsion force simultaneously at the level of the kinetic model. We anticipate 
that similar phase transitions will emerge and will be described by the same limit models, 
with possibly different coefficients. Other points would be worth being developed. The 
computation of diffusive corrections to the SOH model for instance would be of great 
practical use. Further investigations of quantities attached to the social force are also 
very promising, such as the role of a possible anisotropy of the observation kernels, or of 
a different scaling of its range. Finally, the understanding of the transition between the 
two phases and how the two different models can be matched at this transition is also 
crucial for applications. 
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A Appendix: Properties of the collision operator Q 



Proof of Lemma I3.lt Observe: 



M Uf V v (—) 



f J 



A„/ + V„ ■ [-V„ (lnM U/ )/] 



A„/ + V, 



v — u 



T 



f 



Then (13. follows. Formula (I3.19P is a consequence of Green's formula. Now, suppose 
that / G Ker Q, then: 



/ 



dv 



(A.50) 



This implies that there exists a constant C such that / = CM U . In particular, / is of 
the form CM U for a given vector u G M d . Reciprocally, let u G 1R and C a constant, and 
construct / = CM U . Then 



/ fv dv J M u v dv 
J fdv ~ J M u dv 



u. 



Then / = CM Uf . So, by (ETT7jl . Q(f) = and thus / G Ker(Q). This ends the proof of 
Lemma 13.11 ■ 



Proof of Proposition 13.31 We first show that there are d + 1 obvious linearly inde- 
pendent CI's. More precisely, suppose ip{v) = 1 or ip(v) = v. Then ip is a collision 
invariant. Indeed, clearly, such -0 satisfies that ^ G V Uj . for all / such that f/M u . G V^.. 
The statement that ip(v) = 1 is a CI follows from applying (13.181) with g — 1. Let now 
fc G {1, . . . , d} and consider i/j{ v ) = v k- Then 

J Q(f)v k dv = -D j M Uf d Vk (j^-) dv 

= -D J (d v J-fd Vk (\nM Uf ))dv. 

Using Green's Theorem, the first term disappears, and it follows, from the definition of 
Uf, that 



Q{f)v k dv 



f (v k - (u f ) k ) dv = 0, 



which shows that ip{ v ) = v k is also a CI. Thus, C contains a (d + l)-dimensional linear 
space, namely Span{l, v±, . . . , Vd}. We now prove that C is identically equal to this space. 
Let u G M. d . Then, we define the operator R(u; f) as follows: 



R(u; f) = DV V 



M u Vj—) 



(A.51) 



24 



for all / such that f/M u G V u . We notice that for given u G M. d , R(u\ f) is a linear 
operator with respect to / and that 

Q(f) = R(u f ;f). 

Then, ip is a CI of Q if and only if 

J R(uf,f)*l>dv = 0, V/ s.t. f/M Uf G V Uf and tp G K /5 

or equivalently, if and only if: 

VutE* we have: 6 V„ and / Z)^ = 0, 

V/ such that f/M u G 14 and Uf — u. 
As a first step, we fix it G M d and find all ip aV u which satisfy: 

y R(u; f)ipdv = 0, V/ such that f/M u G K and u f = u. (A. 52) 



Then, we will make u arbitrary. We note that any constant ip is a solution of (lA.52j) . 



Saying that Uf = u is equivalent to saying that f f(v — u)dv = 0. So ^ defined by 
(1A.52I) is such that the following implication holds: 

f(v k - u k )dv = 0, VA;G {!,..., 4 [ R(u] f)ij)dv = 0. (A.53) 



Since both sides of the implication of ()A.53j) are linear forms of /, by a standard theorem 
[T3] . there exists /3 U G M d such that 

R(u; f)ipdv = /3 u - / /(u — u)dv. 



Using ( 1A.51I) and Green's formula (see (13.181) ). and introducing the change of function 



g = f /M u , it follows that ip is a solution of the following problem: 

J M u V v gV v iP dv = J g(3' u ■ (v - u) M u dv, Vg G V u , (A.54) 

with j3' u = —j3 u /D. We will drop the primes in the following. 

Let (f) = /3 u -(v — u). Then, problem ( 1A.54|) for ip can be equivalently written according 



to the variational formulation: 

(i>,g)v u = ((i>,g)H u , v<?gK- (a.55) 

Since any constant if; is a solution of ( 1A.52I) . we can subtract J ipM u dv from ip and assume, 
without loss of generality, that ip is such that J ip M u dv = 0. Next, we have: 

(pM u dv = (3 U - J (v - u) M u dv = 0, (A.56) 
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by the definition of M u . Therefore, if g is a constant, we have both (if>,g)y u = (by 
the definition ( I3.16P of (■, -)y u ) and ((f>,g) Hu = (by (1A.56|) ). Therefore, it is possible to 
restrict ( 1A.55I) to functions g such that (g, 1)h u = 0. Consequently, we define the following 
space: 

V u := {(f) eV u : J (f)M u dv = 0}, 
and variational formulation (1A.55jl can be made precise as follows: 

Find iff £ V u such that (if), g) Vu = ((f), g) Hu , V# G V u . (A.57) 
The Poincare inequality for Gaussian measures shows that for all (f> £ V u , 

H 2 K + (0,l)k>Hk- (A.58) 

Then following from (IA.58[) . the bilinear form (•, ■)y u is coercive on V u . By the Lax- 
Milgram theorem, we deduce that there exists a unique ip £ V u such that ( 1A.57I) holds. 
But on the other hand, obvious calculation shows that if)(v) = Tj3 u ■ (v — u) is a particular 
solution. Since it belongs to V u , it is the only solution of the variational formulation 
(007) . 

Adding any constant, we have just shown that any solution of (IA.52|) is of the form 
if)(v) = a + (3 ■ v, where a £ M and £ M. d are arbitrary. Since this form is independent 
of m, we deduce that such -0's are solutions of ( 1A.52I) for all u £ M. d and are therefore the 
only CI's. This ends the proof of Proposition 13.31 ■ 



B Appendix: Expansion of the force term 



Proof of Lemma 15.11 By the change of variables y = x — eR^ and Taylor's formula, 
we have 



K( l -^-^)cuf(y,u),t)dycLu- 



eRf K(\t\)u 



f-eRV x f-i 



(eR) d (j(x,t)+e 2 R 2 kA x j + 0(e 4 )). 



I K(\£\)uf(x -£^,cu,t) d£dcu 
e 2 R 2 



-Dif:(Z®£) + 0(e £ 



x, u), t) d^ du) 

(B.59) 



We have used the definition of k and the evenness of K with respect to £ in order to cancel 
the odd order terms of the expansion. We have denoted by D 2 f the Hessian matrix of 
/ with respect to x (i.e. the matrix of the second order derivatives) and the symbol ':' 
refers to the contracted product of tensors. In a similar way, we have: 



f ~ v\ 

eR 



)f(y t u,t)dydu = (eR) d (p(x,t) + e 2 R 2 k A x p + 0(e 4 )). (B.60) 



Therefore, by expanding the ratio of (1B.59P and f1B.60h up to the fourth order and in view 
of fl2~T0l) . we find f T5~2l) with f l5~26l) . Formula f l5~25l) immediately follows from fl2~T0l) . ■ 
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C Appendix: proof of Theorem 15.2 



Multiplying ( 15.27ft by the collision invariants 1 and v , we are led to the conservation 
equations (in the same way as in Section H]): 

d tP £ + V x -f = 0, (CM) 

d t f + V x • S £ = T~ 1 q £ + eo- x r\ (C.62) 

with p £ , j £ = p £ u e , £ £ , q £ given by (14. 6p and (14.71) and 

r s = j ' u 1 f J £ dv = u 1 fe p £ , (C.63) 
(since u\ e does not depend on v). We also note that we can write 

S £ = p £ (u £ <g> u £ ) + S £ , S £ = [ f £ (v - u £ ) ® (v - u £ ) dv. (C.64) 



The Chap man- Enskog expansion consists in closing the expressions of S £ , q £ and 
r £ by a first order expansion of f £ . To this end, we write the so-called macro- micro 
decomposition: 

f = p £ M u s+efl (C.65) 
By the definition of p £ and u £ , we have 

f £ dv= / p £ M u edv and / f £ vdv= / p £ M u evdv. 



Consequently 

J fl dv = 0, J f{ v dv = 0. (C.66) 

The first term at the right-hand side of (1Q.65D is the macroscopic part, as it is proportional 
to a local thermodynamical equilibrium and carries all information about the moments 
of the solution. The second term is the microscopic part. It carries no information about 
the macroscopic moments but instead carries information about the discrepancy between 
f £ and the local thermodynamical equilibrium. From Section HI we know that, provided 
that p and u satisfy the Euler equations, the microscopic part is small of order e. This is 



why this microscopic part is multiplied by e in ( 1C.65|) . We stress the fact that there is no 
approximation involved (at this step) in ( 1C.65j) : it is a mere definition of f{. 

Inserting (1C.65j) into the formulas providing the expressions of the quantities involved 
in the moment equations (10.61 P . (IC62[) (specifically equations (IC64j) . ( 14 .7p and (1C.63P ). 
we can write the moment equations as follows: 

d tP £ + V, ■ (p £ u £ ) = 0, (C.67) 

d t (p £ u £ ) + V, ■ (p e u £ ® u £ ) + TV x p e = —p £ u £ ( l 1 \ 1 - 1) 

+ e{B\ + B £ 2 + Bl) + 0{e 2 ), (C.68) 
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with 



Bf 



2 = a V e u x pSMue = ^f(2(y x p e ■ V> £ + p £ A x u £ 



B\ = -W x ■ 



fl (v — u e ) ® (v — u e ) dv 



(C.69) 
(C.70) 
(C.71) 



where we have used Green's formula for B\ and ( 15. 26ft . ( 1C.63I) for Bf. Now, in order to 
compute B\ and B$, we need to evaluate ff. But since we only look for 0{e) correction 
terms, and fl is multiplied by e, we may compute fl up to terms of order 0{e). 
Inserting (1C.65P into (I3.17p . the collision operator Q can be written: 



Inserting it in ( 15.27ft . we get: 



sDV, 



(C.72) 



dtf + V, ■ {vD + r- l V v ■ ((1 - ^\)vf) + ea^V v • 



(C.73) 



But we can neglect all terms of order e or more in (1C.73I) . Therefore, we are led to the 
following equation for fl : 



M„.V„(i) 



TV + 0(e). 



where 



TV = d t ( P £ M u e) + V, • {vp E M u e) + T~ 1 V V ■ ((1 - l±-) V p E M u e 



(C.74) 



(C.75) 



The inversion of (10.74j) will give us fl. Once fl is obtained, we insert it into (1O.69j) - (1O.70[) 
and this leads us to the expressions of the Navier-Stokes terms. We first compute the 
right-hand side 7Z £ : 

Lemma C.l. We have (dropping the superscript e for the sake of clarity): 

K = -pM u {h(v - u) : V x u + —[(d+ 2)Tb(v - u) - a 2 dc(v - u) 

+{d + 2)e(v -u):(u®u) + 3T 1/2 (d + 2)g(v -«)•«]} + O(e), (C.76) 



28 



where b(w), c(w) are scalars, e(w), h(w) are tensors and g(w) is a vector, and are given 
by: 

h(w) = Id — - , b(w) = j — ^ (1 - . ,' ' . V c(w) = l-++-, (C.77) 

K i T , v j rp v (d + 2)T J ' v ' dT ' v ' 

e W = ( 1 -(dU)r) W -(dT2)T' ^^"(dT^T^ (C78) 

i?y construction, these quantities are dimensionless. 

Proof. We use the hydrodynamic equations ( 1C.67I) . ( 1C.68I) (dropping the 0(e) terms), 
in order to replace the time derivatives by space derivatives in ( 1C.75I) . This procedure is 
a classical step of any Chapman-Enskog expansion. For simplicity of notation, we omit 
the dependencies on e. Concerning the first two terms of (jC.75p . we write: 

(d t + v- V x )(pM u ) = M u {(3 t + v- V x )p + p(d t + v ■ V x )(lnM u )} 

{v — u 
(dt + v V x )p + p—^- -(d t + v V x )u 

= M u {d t p + u ■ V x p + (v - u) ■ V x p 

v — u 

+P T ■ (d t u + (u ■ V x )u + ((v - u) ■ V x )u) 



Now, using ( 1C.67I) . ( 1C.68I) (dropping the 0(e) terms), we have: 



d t p + u- V x p = -pV x ■ u + 0(e), 

d t u + u- V x u = -TV, lnp - r-^( N2 + ( ^ +2)r - l) + 0(e). 
Inserting these expressions into the previous ones leads to: 

(dt + v- V x )(pM u ) = M u { - pV x -u+(v-u)- V xP 

v-u r _ , , ,\u\ 2 + (d + 2)T \ 
+P— ■ [ ~ TV X In p - r~\{^ L L- - i) 

+((v-u)-V x )u]}+0(e) 
r _ v-u r _j ,\u\ 2 + (d + 2)T . 

u {-pV x . U + p—- [-T ^(J-J K - ^-1) 

+((v-u)-V x )u]} + 0(e) 
u { - Ph(v - u ):V x u-^ ( |M|2 + ^ + 2)T - l)(v -«)•«} + 0(e) 



M 



(C.79) 



29 



For the third term (]C.75p . we compute 



12 



r- l V v ■ ((1 - ^vpMu) 



a 2 



-T 



- VM „{( rf + 2)^ + (l-^p)-rf}. (C.80) 



Collecting the second term at the right-hand side of (1C.79P together with (IC.80P leads to 



TO? 



an expression S which we can split in the following way: 

S = -M u ^- ( H2 + ( ^ + 2)T -l)(v-u).u+ -V, ■ ((1 - ^)vpM u ) 

M u {(d + 2)Tb{v -u)-a 2 dc(v — u) 
+{d + 2)e(v -u):(u®u) + 3T 1/2 (d + 2)g(v - u) ■ u} , (C.81) 
Collecting this expression with the first term at the right-hand side of ( 1C.79I) leads to 

In order to solve equation ( 1C.74I) for /f , we need to solve equations of the type 

L u f := -DV V ■ [M U V„(-^-)] = g, (C.82) 

where u is an arbitrary vector of M d and g is a given function. We refer the reader to 
Section |3~21 for the definitions of the spaces H u , V u and V u . We state the following lemma, 
whose proof is identical to that of proposition 13.31 and is left to the reader. 



Lemma C.2. Let g be such that g/M u e H u . Then, equation ( fC. 8fy) has a solution if 
and only if g satisfies the solvability condition 

gdv = 0. (CM) 

Under this condition, problem ( f(7. 8B) has a unique solution f such that f /M u G V u , or in 
other words, such that f /M u £ V u and satisfies 



fdv = 0. (C.84) 

This unique solution f is denoted by f = L~ Y g and L" 1 is called the pseudo-inverse of 
L u . The set of solutions to equation ( IC 82\) is given by L~ l g + Span(M u ). Additionally, 
if g is such that 

gvdv = 0, (C.85) 



then, f = L^g satisfies 

r fvdv = 0. (C.86) 
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Now, we verify that each of the elementary functions b(v — u)M u , c(v — u)M u , e(v — u)M u , 
h(v — u)M u and g(v — u)M u satisfy both (IC.83P and fl0.85[) and therefore, that the 
corresponding equation ( 1C.82I) is invertible. More precisely, we have: 

Lemma C.3. The functions h{v — u)M u , b(v — u)M u , c(v — u)M u , e(v — u)M u and 
g(v - u)M u satisfy IRTM) and UTM) . We introduce H(v - u)M u = -Z,~ 1 (/i(v - u)M u ) 
and similarly for B, C , E and G. We have: 



~h(w), B{w) = -^{c{w) + - d b{w)), 



H(w)- 

C{w) = ~c(w), E(w) = --e(w), G(w 
and H(v — u)M u through G(v — u)M u satisfy (C.84) and $G86}) . Then: 



fl = a P M u \hi{v-u) :V x u 



TO? L 



(d + 2] 



(d + 2)T u . . dAd + 2)T 2N . . 
e(v - u) : (u (8) u) + T 1/2 (d + 2)g(v - u) ■ u 



(C.87) 
(C.88) 



0(e). (C.89) 



Proof. The proof that h(v — u)M u through g(v — u)M u satisfy (IC.83I) and (IC85[) eas- 



ily follows from classical formulas for moments of the Gaussian, which we leave to the 
reader. Then, we apply Lemma |C. 2 1 which gives the existence of L~ l (h(v — u)M u ) through 
L- l (g(v-u)M u ) and the fact that they satisfy flClMjl and (!Tl86l) . Formulas f|Cl87jl . ( 10881 



follow from explicitly computing the action of L~ l and using the uniqueness statement 
of Lemma IC.2I Finally, equation ( 1C.74|) for ff, which can be written (up to order 0(e) 
terms) —L u f{ = 71, can be solved by ff = —L^IZ since according to the first equation 
of (IC.66p . fl satisfies (IC84[) . By the linearity of L u and the decomposition (1C.76P of 71, 



we can write: 



f{ = -pM u {H(v - u) : V x u + — [(d + 2)TB(v - u) - a 2 dC{v - u) 

+(d + 2)E(v -u):(u®u) + 3T 1/2 (d + 2)G(v -u)-u]\ + O(e). (C.90) 
Thanks to (ITI871) . (!Cl88ll . equation (!Tl89l follows. ■ 



We are now in a position to calculate B\ and B 3 (see (1C69j) . (lC.7ip ). We state: 
Lemma C.4. We have: 

Bi = \p{ [(V, • u)u + V,(^) + [u ■ V x )u] + ^u\^- -v]} + 0(e), (C.91) 

B 3 = /iV x . • (p€(u)) + VMP, u) + XV X ■ (pu ® u) + O(e), (C.92) 
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with £{u), n, A and ir(p,u) are defined at theorem [5H and v is given by: 



Proof. We first consider B\. Splitting v into (v — u) = u and using (1C.66I) . we have 



Bi = — — / {\v - u\ 2 (v - u) + \v - u\ 2 u + 2((v - u) <g> (v - u))u}f e dv 
ra 2 J 

= Ji + J 2 + Js- 

Ji involves an integral of an odd power of v — u and only the g term in ff contributes to 
it. J2 and J3 involve integrals of even powers of v — u and therefore, only the h, b, c and e 
terms need to be taken into account. The computation of these terms rely on computing 
moments of the Gaussian which are left to the reader. We find: 

2T 2 a(d + 2) 
r z a q 

T<r ,„ , daTAd + 2)T , (d + 2)aT . l2 

h = V x -« pM + — P ^ l]pu+- — p\u\\ + 0(e), 

ra 1 r 2 a 2 L a 2 J r^a 4 

Td . , ,rp 2aTAd + 2)T , 

+ 0(e). 



r 2 a 4 



Now, adding up these three expressions leads to ( 1C.91I) . 

We now turn our attention to B3. For this purpose, we compute the tensor 

U = J f{ (v - u £ ) ®{v- u £ ) dv. 

Since the integral involves an even power of v — it, only the h, b, c and e terms of /f need 
to be taken into account. The computation leads to 

U = -^ P ((V x u) + (V x uf) - [^±|^ - l]ld 
2 t a 

-^p(\u\ 2 ld + 2u®u) + 0{e). 
ra 2 



Inserting this expression into ( 1C.71j) leads to ( 1C.92I) . 



Now, by adding the expressions of B\ through B3 found above into ( 1C.68I) . we find the 
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following momentum equation: 



^(pu) + V x • pw (8) u + TV x p = — pu(^ K - '- 1) 



A ( T ,„ \ „ ,\u\ \ . * -| d + 8 Jul 2 n i 

+ e g <°{ K V * ' u ^ u + V * (^") + ( M ' V *H + ~ ^ [^r - ) 

k 

+ e—(2(V xP -V x )u + P A x u) 
a 

+ epN x ■ {p£{u)) + V x 7r(p, u) + AV, • (pu ®u) + 0(e 2 ). (C.94) 

Now, the second term at the left-hand side of (I5.29P combines the second term at the left- 
hand side and the last term of the right-hand side of (10.94j) ; the third term at the left-hand 



side of (15.291) combines the third term at the left-hand side and the penultimate term of 
the right-hand side of (1C.94I) ; the first term at the right-hand side of (15. 29ft combines the 
first term at the right-hand side and the last term of the second line of ( 1C.94I) ; and the 
other terms are unchanged but merely re-ordered. This ends the proof of Theorem 15.21 ■ 
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